Planar Voronoi cells and the failure of Aboav's law 
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Aboav's law is a quantitative expression of the empirical fact that in planar cellular structures 
many-sided cells tend to have few-sided neighbors. This law is nonetheless violated in the most 
widely used model system, viz. the Poisson- Voronoi tessellation. We obtain the correct law for 
this model: Given an n-sided cell, any of its neighbors has on average m n sides where m n = 
4 + 3(n/n)~ 5 + ... in the limit of large n. This expression is quite accurate also for nonasymptotic 
n and we discuss its implications for the analysis of experimental data. 



Two empirical rules play a key role in studies of planar 
cellular systems: Lewis' law and Aboav's law. Both are 
statements about the statistics of a cell's most prominent 
properties, viz. its area and its number of sides. Lewis' 
law says that the average area A n of an n-sided cell 
increases with n as A n = ao(n — no)/X, where ao and uq 
are constants, and A is the two-dimensional cell density. 
Aboav Q noticed that many-sided cells tend to have few- 
sided neighbors and vice versa. He expressed this corre- 
lation in terms of the average m n of the number of sides 
of a cell that neighbors an n-sided cell. Aboav's law, also 
called the Aboav- Weaire [3( law, asserts that 



b 

n 



(1) 



where a and b are constants. This law is in widespread 

use HE00SHE1E1E1E1E3 in the anal y sis ° f 

experimental data on cellular structures. 

In nature, planar cellular systems come in a wide vari- 
ety. They include biological tissues 0, 0> Q , polycrystals 
|2j, cells formed by particles trapped at a water/air inter- 
face , cells in surface-tension driven Benard convection 
0, in two-dimensional soap froths and in magnetic 
liquid froths |9j. 

Alternatively, the cellular structure may appear when 
the data are subjected to the Voronoi construction [Icj . 
Examples are hard disks on an air table a binary 
liquid during late sta ge c oarsening fl2j |. two-dimensional 
colloidal aggregation [l3j, nanostructured cellular layers 
[l4| , and studies of two-dimensional melting 0] . 

Many of these cellular systems, in spite of all their di- 
versity, closely obey Lewis' law from n ss 5 up; also in 
many, Aboav's law appears to hold with good accuracy in 
the full experimentally accessible range, i.e. from n = 3 
to n typically between 9 and 12. Fitting QJ to experi- 
mental data leads to values for the coefficient a ranging 
from 4.6 to 5.3 and for b from 7.0 to 8.5. This observed 
similarity of behavior, sometimes called "universality," is 
generally attributed to the strong geometrical constraints 
that accompany a division of the plane into convex cells; 
by contrast, physical or biological mechanisms are be- 
lieved only to lead to corrections. It is of obvious interest 
to subtract from the experimental data any purely geo- 
metrical effect that one can isolate in order to identify 



the mechanisms at work that cause the corrections. 

Although Lewis' and Aboav's laws arose initially 
merely as good descriptions of the available data, several 
"derivations" have since conferred to them the status of 
theoretical truths valid for all n and embodying a purely 
geometrical theory. Most attempts to derive these 
laws apply maximum entropy principles to a hypothe- 
sized entropy functional ^3 • A critique of this usage of 
the maximum entropy method is due to Chiu |l8j . who 
has shown that no firm conclusions can be drawn from 
it. In any case, none of these derivations connects either 
Lewis' or Aboav's law by a first principle calculation to 
a microscopic model of planar cellular structure. 

First-principle approach. - The simplest microscopic 
model of a planar cell model is the Poisson- Voronoi (PV) 
tessellation: It is obtained by constructing the Voronoi 
cells ^3 °f a configuration of randomly and uniformly 
distributed point centers in the plane, for convenience 
often called "seeds" (but without the implication that 
they are material). The PV tessellation is therefore a 
natural candidate for the purely geometrical theory, and 
certainly the easiest one to handle. 

The statistical properties of the PV tessellation have 
been studied analytically and by Monte Carlo methods 
(see Ref. 0] for a review). Nonetheless, a first princi- 
ple derivation of expressions for A n and m„ has long 
seemed forbiddingly difficult. Numerically it is known 
UlEllii that Aboav's law is not exactly valid for 
the Poisson- Voronoi tessellation: Whereas Eq. pre- 
dicts nm n to be linear in n, this quantity shows in fact a 
very small but distinct downward curvature for the PV 
tessellation. To accommodate this discrepancy modifica- 
tions of Aboav's law have been suggested, in particular, 
by Boots and Murdoch and by LeCaer and Ho [20| . 
The former authors write m n = A-\-Bn~ l +Cn~ 2 , which 
was later found to be true exactly for a class of graphs 
in field theory 1231. 

A recent letter [24( has achieved a step forward in the 
statistical mechanics of planar Voronoi tessellations. It 
opens up the possibility for an expansion in powers of 
n~2 of all quantities of interest related to the n-sided 
Voronoi cell; an immediate result was that Lewis' law 
holds with coefficient ao = \ for asymptotically large n. 
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FIG. 1: Schematic picture of an n-sided Voronoi cell having 
n « 100 and containing a seed at O. First and second neigh- 
bor cells contain seeds at Fi , Fb , F3 , . . . and at Si, S2, S3, . ■ ., 
respectively. All solid line segments separate Voronoi cells. 
Among these, the heavy solid line To is the perimeter of the 
n-sided cell, which is close to circular. The heavy solid line 
Ti separates the first from the second neighbors. Both are 
piecewise linear on a scale n - '. On the scale of order 1 the 
incipient piecewise parabolic structure of Fi is discernible. 



Aboav's law. - Whereas Lewis' law refers to a single 
cell, Aboav's deals with the intrinsically more difficult 
problem of correlations between cells. Here we track 
down the implications of the nT^ expansion for Aboav's 
law. We consider an n-sided cell with n very large. Cells 
with very many sides are extremely rare, but if one oc- 
curs, then its environment must look as depicted in Fig.Q 
where the n-sided cell of a "central" seed at O is sur- 
rounded by n strongly elongated first neighbor cells con- 
taining seeds Fi. Independent evidence for such a geome- 
try comes from work by Lauritsen et al. |22( , who Monte 
Carlo simulated a Hamiltonian favoring the appearance 
of many-sided cells. 

In Fig.QJat least four different length scales play a role, 
each of them proportional to its own power of n. Heavily 
relying on the work of Ref . |2*3 | , which for the present pur- 
pose we extend and interpret, we establish the following 
list of scales: 

(i) The perimeter r of the central cell typically runs 
within an annulus of center O, of radius R c — in/AirX)? , 
and of width of order 1. Hence for n — > 00 the perimeter 
tends towards a circle of radius R c . Consequently in that 
limit the first neighbors Fi will be on a circle of radius 
2R C . 

(ii) Two successive vertices on the perimeter of the 
central cell have an average distance £ VCT t — 27rR c /n = 
(ir/nX)?. Consequently, two successive first neighbor- 
seeds Fi and F i+ i have an average distance 2£ veTt . 

(iii) Locally the positions of the vertices of the cen- 
tral cell are strongly aligned: Their radial coordinates 



have rms deviations of order n~2 with respect to the 
locally averaged radius. For the purpose of the present 
discussion we may set these deviations equal to zero; their 
smallness implies that for n — > 00 the perimeter of the 
central cell becomes a smooth curve. Similar statements 
hold for the curve, not drawn in the Fig.^ that links the 
successive first neighbors. 

(iv) In the region to the right of the heavy solid line 
Ti, which is occupied by second and further neighbors, 
the seed density has its "background" value A, i.e. is of 
order n° H|. 

On the basis of this picture we reason as follows. The 
perimeter of the central cell carries a vertex line density 
Pvert = 1/^vert = (^A/tt)^, which for n — ► oo tends to 
infinity. In spite of this diverging line density the surface 
density A of the seeds to the right of the curve Ti stays of 
order n°. It follows that the central cell will have ~ ni 
second neighbor cells and that each second neighbor Sj 
will be adjacent to ~ -nh first neighbor cells Fi (where we 
call the cells by the names of their seeds). Fig. ^ shows 
that under these circumstances each first neighbor Fi is 
most likely to have itself four neighbors, viz. the central 
n-sided cell, a single second neighbor cell, and two other 
first neighbors, and 

We now focus on the exceptional that have five 
neighbors due to their being adjacent to two second 
neighbors Sj and Sj+i- An example is the cell marked 
Fi in Fig.^ which is adjacent to both S\ and 82- Let 
us denote by the fraction of first neighbors that are 
5-sided. In view of the scaling relations that precede we 
expect that f$ = cn — 3 + • • ■ , where c is a numerical co- 
efficient and the dots indicate terms of higher order in 
n _ 2. Any 6- and higher-sided Fi will contribute only to 
these dot terms. Hence we have 

m n = 4(1 -/ 5 ) + 5/5 

= 4 + cn~3 + . . . . (2) 

We next consider the solid line Ti in Fig.^ which sepa- 
rates the central seed's first from its second neighbors. In 
the large n limit the curve linking successive first neigh- 
bor seeds becomes a circle which may locally be replaced 
with a straight line. In Fig.[3 this same straight line is 
represented by the x axis and the region of space contain- 
ing the second and further neighbors by the half-plane 
y > 0. The second and further neighbor seeds are uni- 
formly distributed in the upper half plane with the back- 
ground density A. Since the first neighbors Fi are dense 
on the x axis, the curve T% now divides the half-plane 
y > into a lower part of points closer to the x axis 
than to any of the seeds, and its complement. Hence 
the function y = Ti(x) is piecewise parabolic; its incipi- 
ent parabolic segments are discernable in Fig.Q] To each 
cusp of Ti(x) corresponds a 5-sided first neighbor cell. 
Let £ cusp be the average distance between the abscissae 
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FIG. 2: The a;-axis represents a continuum of first neighbors. 
The upper half-plane is randomly filled with seeds of uniform 
density A (here A = 1). The piecewise parabolic curve y = 
Ti(x) separates the region of the half- plane closer to the x 
axis than to any of the random seeds from its complement. 

o i 

The abscissae of the cusps of Fi have a density of | A 2 on the 
x axis. Solid and open squares represent second and further 
neighbors, respectively. Note that the horizontal and vertical 
scales are different. 



of two successive cusps of Ti(x). Then it is clear that 
/s = 2^ V ert/^cusp- The determination of £ CU sp for given 
seed density A in the upper half plane is a problem in 
statistics that yields £ CU sp = 2/(3 As) [2(|. Using the 
expression for £ VC rt found above we therefore have that 
fs = 2(7r/?iA)2 x 3A3/2 = 3(71/71)2, whence c = 37T2. 
Substituting this in J2J we conclude that for the Poisson- 
Voronoi tessellation m n is exactly given by 



m n = 4 + Z\ — 
n 



(3) 



The inverse square root decay and the limiting value 
TOoo = 4 of Eq. J2J are in contradistinction to Aboav's 
law Q. This equation explains for the first time the 
downward curvature observed |l9ll20ll2lLl22l | in the nm n 
vs. n curves for the PV tessellation. 

Data analysis. - We cannot assess a priori the applica- 
bility of Eq. © to m n data in the range of n covered by 
experiments and simulations. We have plotted in a new 
way in Fig. [3] the Monte Carlo data for m n due to Boots 
and Murdoch • The dotted curve shows the best two- 
parameter fit to the data provided by Aboav's law 0J; 
it is obtained for a = 5.251 and b = 5.755 [lfjj and is 
generally taken as evidence that this law fails for the PV 
tessellation. The dashed line in Fig.[3]represents the right 
hand side of Eq. @ with the dot terms neglected. It ap- 
pears that in the regime of the data our zero-parameter 
asymptotic result is numerically only slightly more off 
than Aboav's two-parameter law [2]]]. This high degree 
of accuracy is remarkable in view of the fact that Eq. © 



FIG. 3: The two-cell correlation function of n. Open 

squares: Simulation values due to Boots and Murdoch |l9j|. 
Dotted line: Best two-parameter fit with Aboav's law Q. 
Dashed line: Our zero-parameter exact asymptotic result (J^J , 
which in this representation is a constant equal to . 



rests on an expansion around the improbable event of 
a large n-sided cell. But more importantly, there is full 
compatibility between our dashed line and the downward 
trend of the data for n larger than about 10. One may 
attempt to improve our Eq. by including higher order 
terms in n~2 in the series, the simplest possibility being 



m. 



In d 
•3 A /- + - 

n n 



(4) 



where d is adjustable. However, there is no reason to 
expect that a truncated expression, whether Q, |QJ, or 
yet another one, will yield the exact m n for all n. 

The PV tessellation serves as a model of reference for 
general cellular systems much in the same way as an ideal 
gas does for real gases. The inverse square root decay 
of its two-cell correlation rules out Aboav's law and 
may be viewed as defining a universality class. We can- 
not exclude that someday some other microscopic model 
will be shown to asymptotically obey Aboav's law (JTJ to 
leading order in n~ l . Such a model would then define a 
second universality class. It is unknown at present which 
experimental systems are in the PV universality class. 

Many different factors may potentially lead to depar- 
tures from PV-like behavior, if not to a change of uni- 
versality class. E.g., (i) in some cellular structures the 
"seeds" represent actual particles or larger physical en- 
tities with mutual interactions; (ii) a general planar cell 
structure cannot be derived from a set of point centers 
by means of the Voronoi construction; and (iii) some sys- 
tems, like soap froths, are not in equilibrium but rather 
in a - supposedly scale invariant - coarsening state. One 
may now hope for such cases to be able to discuss the 
deviations from PV tessellation statistics perturbatively. 
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To show how such a discussion might proceed, we con- 
sider briefly and heuristically an example from class (ii) 
above, viz. the Voronoi tessellations associated with hard 
core particles of finite diameter a. The preceding analy- 
sis of the large n limit goes through as long as n < n* , 
where the crossover value n* is determined by the con- 
dition that the distance between adjacent first neighbors 
become comparable to the particle diameter. This gives 
2ttR c ~ n*a, whence n* ~ 7r/(Aa 2 ). For n^n* the repul- 
sion between the particles combined with the condition 
that they be locally aligned imposes that the radius of 
the central cell grows as R c ~ na and hence that /g 
must saturate at a value / 5 = cqoXz where Co is an un- 
known numerical coefficient. Upon assuming that this 
knowledge may be expressed with the aid of a scaling 
function M. we obtain for hard core particles of diame- 
ter a the relation m n = 4 + aX? Ai(nXa 2 ), valid in the 
scaling limit n — -> oo, Xa 2 — > with x — nXa 2 fixed, and 
where the scaling function M. satisfies Ai(0) = cq and 
M.{x) ~ 3(ir/x)z for x — > oo. Hence for hard core parti- 
cles the limiting value of 10 is changed, but the 
decay law remains. Future work will have to deal with 
this and other instances of deviations from PV statistics. 

Conclusion. - The simplest and most widely used 
model of a cellular structure - the Poisson- Voronoi tes- 
sellation - does not obey the most widely used law - 
Aboav's law 0) - governing the two-cell correlations. We 
have shown that for the PV tessellation the n _1 decay 
of Aboav's law should in fact be replaced with an 
decay. Experimental data of sufficient precision and/or 
covering a large enough range should be able to distin- 
guish between the two and thereby shed light on the ques- 
tion of the universality classes and of the underlying cell 
formation mechanisms. We therefore advocate that ex- 
perimental results be analyzed not only by Aboav's two- 
parameter fit l|T)l. but also by the simpler one-parameter 
formula Q , of which the first two terms have a firm the- 
oretical basis. 
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